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ABSTRACT 


The general problem presented in this paper 

is one of estimating the state vector x from the state 

• ■ * 

equation h=Ax where h, A, and x are all stochastic.. 
Specifically, the problem is for an autonomous Martian 
Roving Vehicle to utilize laser measurements in estimating 

i 

the gradient of the terrain. Error exists due to two 
factors - surface roughness and instrumental measurements. 
The errors in slope depend on the standard deviations of 
these noise factors. Numerically, the error in gradient 
is expressed as a function of instrumental inaccuracies.. 
Certain guidlines for the accuracy of permissable gradient 
must be set.. It is found that present technology can 
meet these guidlines.. 



PART 1 


. > INTRODUCTION 

■ ‘ X:- ; r 

A comprehensive navigation system will be needed 

for a proposed Mars rover to safely traverse the surface 

of Mars with reasonable speed, over a long distance.. 

Because the roundtrip communication time to earth requires 

more than 40 minutes, the vehicle’s terrain modeling and 

path selection systems must be autonomous.: The system 

is designed to collect terrain data within a 3 to 30 meter 

range. The range finder, which locates a point on the 

terrain, is a laser/dectector which gives a range 

measurement R^ azimuth angle 0, and elevation angle J5* 

Two- points along the path of the vehicle determine an 

in-path slope, while those across the path compute a cross 

path slope. However, inaccuracy in measurement can 

introduce very large errors . in the computed slopes and 

heights,, which are the main factors in path selections*. 

There are some threshold values for these factors above 

2 

which a change of path is required. 


1 



PART 2 


METHOD OF APPROACH 


From the measurement data it is desirable to 
obtain the maximum slopes and the elevations of the terrain 
in front of the vehicle*, . • ■ 


A*. Transformation of Coordinate Systems . I 

The quantities R, 0 , and J*> are measured with 
respect to the coordinate system h", a", and b M , fixed to 
the vehicle* (see Fig. 1) With laser height at 3 meters 


we have: 


h" r = 3 - Rsin0 
a" = Rcos/}sin8 
b" = Rcos^3cos0 


(la) 

(lb ) 

(lc ) 


The body-bound axis rolls with the angle <j) and 

pitches through an angle f about a reference frame 

h, a, and b formed by the local vertical and' an axis in a 

plane containing the heading and the local vertical. The 

. 3 . 

coordinate transformation is: *, 


(2a) 


h 


h" 

a 

b 

= C(4») B ( ^ ) 

a" 

b" 


where 


C(<|>) = 


cosc^j ,-sin$ 0 
sin<f> cos<|> 0 
0 0 1 


(2b) 





BCO = 


w 


cos? 0 

sin? 

0 1 

0 

-sinf* 0 

coss? 


(2c ) 


B*. Determination of Slopes and Gradient 

A number of measurement points within a certain 
small area of surface, say 0.5m by 0.5m, can.be used to 
determine a plane in space, which may be written as 

h = ax^ + bx£ + x 3 (3 ) 

where x^ and are two constant parameters to be 
determined. Taking the differential of equation (3) we get 

dh = ^da + |rdb = x.da + x 0 db 
da db 1 2 


where, 3h 

x. = — cross-path slope 

x aa 

dh . , . n 4,5 

X 2 = "ko ln “P a ^^ s l°P e 

The corresponding gradient of the plane, which is defined 
as the gradient of the terrain in that small region On the 
planet's surface, is 

? ? i/ 2 

Sg = (x^ + *2 ) . , (4) 

If this slope is less than a predetermined criterion, it 
is considered to be safe for the vehicle to travel ahead. 

In order to locate the plane by a number of 
measurement points, one may rewrite Eq. (3) as 

h i t a i x l + b i X 2 + X 3 <5> 

where h^, a^, and b^ are found from , 0^, /3^ in equations 

th 

(1) and (2 ) for each i point. Theoretically, three 



points determine a plane (i = 1,2,3). For greater \ 
accuracy however,, more than three measurements (i = 1,2,3 
(probably n = 4 or 6) are needed to determine the slopes. 
A complete picture. of the terrain in front of the vehicle 
can be constructed by modeling numerous adjacent planep, 
each covering a small area of surface. 



PART 3 


ANALYTIC SOLUTIONS FOR PARAMETER ESTIMATION 

.In this section we list the solutions of the least 
square estimate, the covariance matrices, and the minimum 
variance estimate when n>3. 


A.. Least Square Estimate of Slopes ; 

If a ^= a j L an< 3 b^=b^ are assumed to be true in 

Eq.. (5)', a least square error estimate can be performed 

which minimizes - 

n * _ o 

( 6 ) 


n 

z 

i=l 


(h.-hp 2 


where h^ is the actual measured height and h^ is the 
corresponding height in the modeled plane.. 

In matrix-vector notation, Eq„ (5) is 

C h = Ax (17 ) 

where 


h = (h, ,h„ , . . . ,h ) T 
1 1 n 


x = (x^ ,x 2 ,x^ ) 


and 


A = 


a l b l 1 
a 2 b 2 1 


a b 1 
n n 


n = 4 or 6 


The least square estimate of the parameter x becomes 

X = (A T A)“ 1 A T h . 


( 8 ) 


B. Perturbation of the Variables 

The least square estimate in the previous section 

7 

assumes the matrix A to be completely deterministic • 


1 6 


7 


In reality, however, there is error involved in the 
determination of a, b, and h, due to the errors in our 
measured values of 4>, f , R, f3 , and 0 . 

If the symbol 8 denotes a perturbation, then 
fih, fra, and $b in terms of £$>,. $$ , £R, Sj3 , SB , are 


fih* 

fia 

= D(h",a",b", 4> } 

S 4> 

+ C(«)B(^)G(R,e,3) 

SR 

S/3 

fib 



_se 


(9a) 


where 

D(h n , a" ,b" ,<}>,^ ) . 

C-h M sin<j>cos^ -a ,f cos4 > “b"sin<j>sin^ ) (-h ,, cos<j)sin^+b ,, cos4>cos^. ) 
( h"cos<J>cos^ -a ,, sin^>+b"cos<j>sini? ) (-h n sin^sin^+b u sin4>cosf ) 
0 


( : -h"cos? -b"sin^ ) 


and 


G(H ,&,#)■ = 


(-sxnfo ) 
(cos;3sin0 ) 
|_(cos^8cos0 ) 


(~Rcos^) 0 

C-Rsip£sin& ) ( Rcos$cos&) 

(-RsinficosB) C-Rcos/3sin0 ) 


(9c) 


The derivation of these matrices are given in Appendix (A)«- 
C.. Covariance Matrix of the Variables 


variables as 


We can define the covariance matrix of the 
8 


( 

fih 

[fih fia fib] 

M = E( 

fia: 

. ■ ' 

( 

fib 

j 


( 10 ) 


8 


where E denotes expected value. 

If 54>, , £R, Sfi , and SB are not correlated^ 

then from Appendix B we have 


(9b) 



M = D 


8 


e ($ 4>) 2 o 


E(S R ) 2 0 0 

0 E(8^) 2 

D T + CBG 

6 E(S/i) 2 0 

o 



o - o e ( 8 e ) 


g t b t c t 


Prom Eqs . (10) and (11) the standard deviations of h, a, 
and b, can be computed in terms of those of <f>, , R, /3 , 

and 0 for each point.; These are known quantities that' 
depend upon the accuracy of the measuring devices, 


, .1 


D. Covariance Matrix of the Slopes 
Eq. (5) can be written as 


( 11 ) 



h = Ax 


(12a) 

If we set h 

= h + 6 h , A = A + 

5 A, and x = x + gx, 

i 

(12b) 

then with the 

aid of Eq. (7) we 

get 



$h = Agx + SAx 


(12c) 

The estimate 

is^ gx = F($h - 

SAx) 

(13) 

where 

„ T-Ty-l-T-T 

F = (A A) A 


(14) 


The covariance matrix of the slopes is determined in 
Appendix C as 

e[SxSx^[ = F[E[sh&h T ] -E [SAxSt 11 ] -E [fih ( SAx ^ +E [gAx ( SAx ) T JJ F T (15) 
where 


A 

Sx = 

S*! 

s* 2 

• Sh = 

Sh l 

Ax = 

^ a l x l‘ , 'k b l x 2 ^ 

• 

# 

* 


/\ 

SX 3 

- 


8h 

n 


(5a n x 1+ Jb n x 2 ) 


Since A is in terms of a and b, then $h and £A can be 
expressed as functions of > $f ? 8 R ? S3? and 80 as given 
in Eq. (9a). Eq. (15) can be evaluated as shown in 
Appendix D. 



E. Variance of the Gradient 


If the symbol °*f denotes the standard deviation 


of Sg, then from Eq. (4) we have 

dSg = (x^+x^ ) -ly/ ^x^dx 1 + (x^+x^ )*" 2/ ^x 2 dx 2 

9 

The variance of the gradient is 


rr 1 - X 1 ^ X 2 7 . ~ X 1 X 2 Z 

°S 3 ~ -2 ~2 + -2 -2 O' + 2 _2 -2 c x | x 1 

J x i +x 2 x l +x 2 2 v i x 


:i7) 


x-^+x 2 


where CT^ = EtSx^) 2 , O^ 2 = S(Sx 2 ) 2 , C^= ECS^gx^) 

2 9 • 2 

The covariances CT V , QT V » and CT V 7) can be found from 

A| *2. P'l^Z. 

Eq. (15). ' , 

The value of gives a rough estimate of the 

accuracy in our estimation of the gradient. If the 

C 6 

estimated Sg is 20 with 0^=2 , then we can be 68% sure 
that the actual slope is between 18° and 22° and 95% sure^ 
that it is between 16° and 24°. This is important when 
the estimation is close to the. maximum permissable value. 


For the case where x^ = x 2 one obtains from 


Eq. (4) 
Thus 


Sg = /2 x 1 

dSg - \fz d'A 1 and Q^= 2C* ( 

8 

Minimum Variance Estimate 


(17a) 


■ A generalized form of Eq. (7). is 

v h - Ax (18 ) 

where v is the noise due to measurement of h and A in 
Eq. (5 ) . 

If the expected values of Eq. (18) are 



then 


IQ 

E(h) = h , E(a) = A '(19) 

E(v) = v = h - Ax = 0 (20) 

• l 

by virtue of Eq.. (12c). Subtracting Eq. (20) from 
Eq. (18) with the aid of the definitions of Eq* (12b), . 
one obtains ' . ' 

v S $h - $Ax (21) 

Thus the covariance of v becomes . 

R = Efvv^rr E[(Sh — S Ax ) ( £ h - S Ax )^J 

= E [ShSh^ - E^AxSh^ - E[sh.(5Ax): T ] + e[sAx(5Ax) T ] 

which is the same as those terms inside the bracket in 
Eq* .(15)* These terms are evaluated /in Appendix D* 

The cost functional J can be expressed as 

J - v^R*~^v s (h~Ax)' T R ^(h-Ax) 

T —1 _rn T — 1 T —1 — — T T — 1 — 

= h R h-x A R h-h R Ax+x A R Ax (23) 

The minimum variance estimate of x can be obtained by 

taking the minimum of J with respect to x. 

Min- J = -r=r = -2A T R~ 1 h + 2A T R _1 Ax = 0 (24) 

x dx , 

Thus the minimum variance estimate x becomes 

x = (A T R“ 1 A): _1 A T R“ 1 h (25) 

where R can be obtained from Eq. (22). 

The minimum variance estimate given by Eq. (25) 

is a weighted least square estimate. It will provide a 

better estimate of the gradient then the standard least 

square estimate of Eq. (8) because it gives correct 

9 

weight to the individual measurement points. The minimum 


( 22 ) 



variance estimate utilizes the value x given by the least 
square estimate to evaluate the covariance matrix R given 
by Eq. (22). 

However, the numerical results in the example 
computed from Eq. (25)‘ is very close to the results given 
by Eq. (8) because the weighting factor in the covariance 
matrix R is nearly proportional to an identity matrix. 

That is to say, the weighting factors in the diagonal 
terms of R are equal. This is equivalent to an unweighted 
least square estimate. The details of the derivation are 
shown in Appendix E where R=kl. 

G. Determination of Gradient 

Once the values for the cross-path and in-path 
slopes are determined from the minimum variance estimate, 
the gradient can be determined from Eq. (4). Knowing the 
value of’ the gradient and its variance, the rover can make 
decisions concerning safe and unsafe terrain. 


PART 4 

NUMERICAL RESULTS 

It is assumed that the rover will use a split 
beam where A is the difference in the elevation angle 
between the two beams and A© is the difference in 
azimuth angle between anY two laser pulses. This is 
shown in Fig. 2. 

The values of @ and 0 for each data point, along 
with the magnitude of the cross-path and in-path slopes of 
the terrain, determine the data point spacing^ and the 
distance of the data points from the vehicle. After a 
number of data points are measured, they are transformed 
to the non-rotating reference system by Eqs. (la to 2c). 
Then a least square estimate of the slopes is performed 
which is given by Eq„ (8). This least square estimate 
depends directly upon the values of h,a, and b for each 

• A 

measurement point. Utilizing the value of the vector x 
found in Eq. (8), the values of (j), ^ , R, ,0 , and 6 for 
each data point, and the standard deviations of . these 
quantities, we can ultimately find the magnitude of the 
gradient and its variance thru the proceedures outlined 
in sections 3b-3f. 

The variance of the gradient depends in part 
upon the magnitude of the cross-path and in-path slopes 
of the terrain and the accuracy of the measuring devices. 
It also depends upon ^ , f , H,^3 , and 0 , or alternately, 




upon the roll and pitch angles, the data point spacing, and 
the distance of the data points from the vehicle. How 
much each of these factors affects the standard deviation 
of the gradient is to be found. 

The standard deviation in gradient for any 
combination of these con be determined by utilizing 
Eqs. (9b)., (9c). , (11 ) , (14), (15 ) ,_ and (17). The data point 
spacing and the distance, of the data points from the 
vehicle can be varied by changing the. values of and 6 for 
each data point. The vector x completely determines the 
gradient of the plane we are •measuring’ and along with, 
specified values of 4* , f , Ji> , and 0 we can geometrically 
find the values for R, h,a, b, h"', a", and b" which are 
utilized in above equations. 

- For simplicity , we first set ^ = 0** for all 

the measurement points and assume. the reference gradient 
of the measured plane to be 0° wherever it may be located. 

By applying Eqs. (9b). (9c), (11), (14 X, (15), and (17a) 

and setting 0^ = CT & -■ „l f 5 and G^> 1° in Eq« H we find 

4 

that for any reasonable data point spacing at 20 - 30 meters 

Uk O 

from the vehicle, Oc, ". as on the order of 30 ~ 60 . On 

the other hand, if CQ - <X^ = 0“ in Eq„ (11),. decreased 

to only 2° -3° at that distance. • 

. Since it can be expected that 0$ and are really 
in the neighborhood of 1 ° (due to the constant rock and 
roll of the vehicle as it traverses the surface), the rover 



will have to be provided with a> ’rapid scan* laser.. If 
the scan rate is on. the order of milliseconds, which is 
perfectly feasible with electronic scanning]: then each 
set of 4 adjacent data points is measured practically 
..instantaneously, since the rover motion is on the order of 
seconds. Each of the 4 data points will have. the same 
value for <p and ^ and therefore all 4 points will retain 
the same relative position to each other when they are 
transformed from the vehicle, coordinate system to the 
fixed system. Consequently, the rover can model the 
planes in the vehicle coordinate system and then transform 
. the planes to the non-rotating system, instead of first 
transforming each point to the non-rotating frame. In 
the analysis, this corresponds to setting ^ = 0°in 

Eqs. (9b), (9c), and (11 )' and setting h"=h, a"=a, and 

b"=b in these equations. Since the maximum positive or 

. o 

negative slope that the vehicle can navigate isl25 , 

I . o 

relative slopes might be as high as —50 and still be 
navigable. This is shown in Fig. 3. Therefore,, in my 

O 

analysis, I consider a range of slope changes from +50 
to -50° where possible. 

Also, due to rapid scan, the error introduced 
by a standard deviation in Cfyand c will be consistently 
between 1° and 2° because it only involves the error in 
■■ transforming the already modeled plane from the. vehicle . 
system of reference to the fixed system. Consequently, 
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Fig. 3 Relative in-Path Slopes 



the effect of the standard deviations in <j> and ^ is not 
considered in the analysis and 0 $ and 0 ^ are set to 0 
in Eq.(ll). It must also be realized that any. relative 
slopes between +50° and —50° can represent critical values 
.in the fixed system and therefore all relative slopes are 
important.. 

If the vector x=0 in Sqs.-(15) and (17a), 0 * 0 = 1 ’ 

and <T R =5cm in Eq.(ll), and A /3 and AS (and consequently 
the data point spacing! are varied for each set of 4 points, 
then the graph shown in Fig. 4 results* This corresponds 
tp collecting data from a flat terrain. Each solid line 
represents constant values for A/3 and A0. The quantities 
Ab and/ia are the data point spacings along the in-path 
and cross-path directions respectively. By looking at 
any solid line it can be seen that for any constant values 
of A# and A0 , the quantities Ab and Aa decrease very 
rapidly as we scan closer to the vehicle and consequently 
Ogg rises very rapidly. 

In' choosing an ’optimum’ spacing it must be 

I • • 

kept in mind that by increasing the spacing <\b and Aa, 
the standard deviation 0*53 decreases. This also renders' 
the data less meaningful, as more terrain has been overlooked 
At the 30 meter range, a general picture of the terrain 
with a spacing of 2-3 meters Is sufficient. At close 
range, (4-7 meters) the points should be at least as close 
as 0.66 meters because this is the width of the widest 





2 

navigable crevice. In Fig. 4, the dotted line represents 
CTgg vs. distance from the vehicle for such a scheme, 
where the data point spacing varies from about .5 meters 
at a distance of 4 meters from the vehicle, to a spacing 
of about 3 meters at a distance of 30 meters from the 
vehicle. . 

Each of the plots in Figs .5 and 6 utilize the 
same optimum scheme as the dotted line in Fig. 4. Once 
again the vector x=0 (corresponding to flat terrain).. In 
Fig. 5, is kept at 5cm inEq.(ll) and 0^ and <y & are set 
to 0 , .017 ( 1 * ) , and .1 °( 6 * ) „ It is seen that ^ is 
markedly reduced by decreasing Op, and (Tq from 6' to 1 ' , 
but any further increases in accuracy w ill result in 
bulkier and heavier equipment which is not very beneficial 
By comparing the plots for C^=Oq= 1' and 0 %= cr e =0°, it 
is seen that Oq and Op have very little effect on 0$ 

NJ 

when they have a value of 1’. In Fig. 6, CT R is varied from 
icm to 10cm in Eq.(ll). while CT^ and Oe are kept constant . 
at l'C.,017 ). The quantity (^increases rapidly for each 
value of CT R at close range in Fig. 6. . 

Fig. 7 plots o^ VSo relative in-path slopes of 
from —30° to +50° at 4 meters from the vehicle. This 
means that the vector x is varied in Eqc-(15) and (17) such 
that the cross path slope remains 0? the in-path slope 

— 0 o — , 

x 2 varies from -30 to +50 , and x^ changes m such a 
manner that there is always at least 1 data point with a 
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height h of 0 meters in each group of 4 points. The graph 

• ' o 

only extends to -30 because the laser beam cannot ’see' 
over a negative in-path slope greater than this value at 
4 meters distance. In this graph , Cp = 0^ =1 1 in Eq.(ll) 
and is varied from 1cm to 5cm. The data' point spacing 
is .6 meters.. Fig. 8 is also a graph of 0£^vs. relative 
in-path slqpe from -6° to +50° at 20 meters from the vehicle. 
All the other conditions are the same as in the last graph 
except that the data point spacing here is 1.2 meters 
instead of 0.6 meters. Finally, Fig. 9 plots graphs, of 

O O 

vs. relative cross-path slopes of from 0 to +50 and for 
distances of 4 meters and 20 meters' from the vehicle. Once 
again Op = =1’ and is varied from 2 to 10 cm. 

The data point spacing at 20 meters , from the vehicle is 
1.2 meters and the spacing at 4 meters from the vehicle 
is 0.6 meters. 

Present technology cannot improve CTr below 1cm, 
and even ohis figure is quite lowl^ It is obvious from 
Fig. 7 that Q r R should be as close to 1 cm as possible in 
order that °y will be within an acceptable range at 4 
meters distance and high values of relative in-path slopes. 
Cross-path slopes do not present as much of a- problem and 
a CT R of 2 cm is quite acceptable for the worst cases. 




Standard Deviation in Gradient 
CJ" s g (degrees) 


At 4 m Distance • 

Or= 0.03 and 0.02 m 
Data Point Spacing 0.6m 

At 20m Distance ! 

0" R = 0. 1, 0.05,0.03, and 0.02 m 
Data Point Spacing 1.2 m . 


4m; (7r= 0.03m 


4m; 0 "r= 0.02m 


20m;0'R = 0.lm 

^20m ;(T r = 0.05 m 
-20m;(T R = 0.03m 
‘20m;0 " r = 0.02 m 


Relative Cross -Path Slope (degrees) 
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PART 5 


CONCLUSION 

. "V • . ' • ' ' 

\ A 

This paper suggests a solution to the general 
problem of estimating the state vector x from the equation 
h=Ax where the measurements h and A are all stochastic. 

The proceedures outlined here are applied to a specific 
case in estimating the gradient of the terrain by laser 
measurements from an autonomous Martian Roving Vehicle. 

Due to terrain irregularity and instrumental uncertainties, 
the four measurement points do not fall in the same plane.. 
First, a least square estimate is. performed which assumes 
the matrix A to be deterministic. A. minimum variance 
estimate, which takes into account the error in the matrix 
A as well as h, gives approximately the same numerical 
results because the weighting factors are nearly equal 
to one ancther. The error involved in these algorithms 
can also 1 e estimated and a complete error analysis has 
been presented.. 

# 

| The variance of the gradient should be as low as 

possible but an upper bound of 2 degrees has been set. 

Figs. 4-9 show that the standard deviation of the 

gradient is dependent upon the distance of the modeled plan 
away from the vehicle, the data point spacing, the gradient 
of the terrain, and the values for CX R , , and <y Q . By 

using a rapid scan laser, we do not have to take into 
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consideration the error due to the motion of the vehicle. 
If Oft and are reduced below 1 arc minute, bulkier and 

heavier equipment would result.. In fact, with0£=Oe of 
1 arc minute very little error is introduced into the- 
calculations. Most of the errors come from . In 

order to keep the upper bound of within 2 degrees, the 

quantity 0"^ must be as low as 1 cm. 
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APPENDIX A 

DERIVATION OF THE PERTURBATIONS 


To prove the result in Eq.. (9), the coordinates 
in Eq,. (2) are first perturbed as 


a, 

6h 


"h" * 


‘ h" ’ 


fih" * 

.. fi,a- 

6C(^)B(?) 

a" 

+ c( 4) bU) 

a" 

+C(4»)B(C) 

6 a" 
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b" _ 


. B" . 


fib" _ 


( A-l) 


where 
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and 
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(-cos£fi£) 

0 
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(-cos 5) 

0 ( 

-sinOj 


4,-nl . 


(A-2) 




(A-3) 


perturbing Eq. (la) - (lc) 


lh"' 


’8 r" 

6 a" 

= • 

g(r,o,jS) 

S/3 

fib' 


Se 
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*■ 


(-sin B) 



(-Rcos/3) 0 f 

(cos/5sin8) (-RsiniSsinO ) ( Rcos/3’cos 9 ) 

(cos0cos8) (-Rsin0cos0 ) C-Rcos/3sin0 ) 

By substituting (A-4),. (A-3) and (A-2) into (A-l) one 
obtains 



SR 


S/3 


60 


(A-4) 


Sh" 


r 

h" 

• 

h" 

-) 

r ” 


SR 

fia 

P- 


c 1 («J> ) B( ^ ) 

a" 

b" 

*c ( 4>) b 1 ( ^ ) 

• 

a" 

b" 

j 

S4> 

Sf 

md 

+C(4>)B(f)G(R,&,e) 

t * r \ 

S/3 

JBG^ 



29 


Compart the above Eq. with Eq. (9a), to see that 
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and' 
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0 .. 

The second terra of (A-6) is determined from (2b) and(A-3) 
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By substituting Eqs. (A— 7) and (A-8) into (A— 6) one obtains 
the results in Eq. (9b).. 



APPENDIX B 


DERIVATION OF THE COVARIANCE MATRICES 
OF THE VARIABLES 


To prove the result in Eq. (11) we multiply Eq.(9) 
by its transpose.. 
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Taking the expected value of (B-l) one obtains 
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Since $4 , Sf , SR, S/3 , S3 are not correlated, ( B-2 ) 
becomes Eq. (11) in the text. 



APPENDIX C 


DERIVATION OF SLOPE COVARIANCE (I) 

Here we will prove Eq. (14) and Eq. (15). SJ 
Substituting h=h+£h, A=A+SA, and x=x+5x into Eq.(12l 
we have 

h+Sh= (A+5A) (x+6x:)=Ax + £Ax + A5x + SA Sx (C-l) 

Since h=Ax, and neglecting the second order term 6A8x, 

Eq. (C-l ) becomes 

A Sx = Sh-5Ax . 

— T 

by premultiplying the above equation by A 

A T Agx = A T (Sh-SAx). 

Thus the estimate of 6x becomes 

Sx = (A T A)" 1 A T (Sh-SAx) = F(Sh-SAx) ( C— 2 ) 
which is equivalent to Eqs.(13) and (14). Multiply 
Eq. (C-2) by its transpose to obtain 
6xgx T = F(Sh-SAx) (6h-$Ax) T F T 

= F[5hSh T -<5^h T -Sh(SAx) T +SAx(SAx) T ]F r -. 

( C— 3 ) 

The expected value of Eq. (C-3 ) becomes Eq. (150 in the text 



APPENDIX D 


DERIVATION OF SLOPE COVARIANCES (Z) 

Replacing, term by term, Eq. (16) into Eq. (15) 


we get 
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The different measurement points are uncorrelated. Therefore > the 
off diagonal terms are equal to zero. Eq. (D-2) reduces to 
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Similarly: 
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APPENDIX E 

THE NOISE COVARIANCE MATRIX ' 

From the quantities of gh and gAx given in Eq.(16) 
one obtains for Eq„ (21) ... 

Sh ]L -Sa 1 

v = : 

* 

J h 4 -8 a 4 

Since the measurement points are close together, the 
expected values of the errors are nearly the same, i.e„, 
Etfhj.') = S(8h 2 ) = E(5h 3 ) =? E(8h 4 ) 

EUa^ * E(8a 2 ) ~ E(ga 3 > S E(ga 4 ) 

ECgbj^) ST E(gb 2 ) * E ( 8i> 3 ) = E ( 5 b 4 ) (E-2) 

similarly 

E(£h 2 ) 2 =r E(5h 3 ) 2 ^ E(gh 4 ) 2 
^(ga^ 2 * E(ga 2 ) 2 5 E(ga 3 ) 2 * E(ga 4 ) 2 

ECS^ ) 2 s E(5b 2 ) 2 = E(5b 3 ) 2 =r E(gb 4 ) 2 (E-3) 

From Eq«> (D-3) we have 

E^Sh^] = E(5h 1 ) 2 I ( E-4 ) 

Similarly from Eqs,. (D-4) and (D-5) one obtains". 

E[gAxSh T ] = E [gh ( S.Ax ) T ]J 

«[E(6a Sh^)^ + E(Sb 1 8h 1 )x 2 ] I ( E-5 ) 
Moreover, Eq*. (D-6) results 

e[$Ax(SAx) T ] =[E(8a 1 ) 2 x 2 + 2E(8a 1 5b 1 )x 1 x 2 + Etgb^x 2 ]! 

(E— 6 ) 

Substituting Eqs. (E-4), (E-5) and (E-6) into Eq. (22) gives 



R = kl 


( E— 7 ) 
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where k = E($h^ )^-2E(Sa^$h^ )x^-2E(6b^Sh^ )x^ 

+ E(Sa 1 )^x^ + 2E(5a^§b^ Jx^x^ + ECSb^ (E-8.) 

which is the predicted result in section 3f. 



